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Abstract 

We discuss same-sign dilepton production mediated by Majorana neutrinos in high- 
energy proton-proton collisions pp i^i'^X for i, i' = e, fi, t at the LHC energy 
y/s = 14 TeV, and in the rare decays of i^", D, Dg, and B mesons of the type M"*" 
M'~i~^i'~^. For the pp reaction, assuming one heavy Majorana neutrino of mass rriN, 
we present discovery limits in the {rriN, lUmUi'Nl) plane where UiN are the mixing 
parameters. Taking into account the present limits from low energy experiments, 
we show that at LHC for the nominal luminosity L = 100 fb~^ there is no room for 
observable same-sign dilepton signals. However, increasing the integrated luminosity 
by a factor 30, one will have sensitivity to heavy Majorana neutrinos up to a mass 
rriN ^ 1.5 TeV only in the dilepton channels fifi and fir, but other dilepton states 
will not be detectable due to the already existing strong constraints. We work out a 
large number of rare meson decays, both for the light and heavy Majorana neutrino 
scenarios, and argue that the present experimental bounds on the branching ratios 
are too weak to set reasonable limits on the effective Majorana masses. 



1 Introduction 



Recent results from the KEK to Kamioka long baseline neutrino experiment (K2K) 
jl] strengthen the neutrino oscillation interpretation of the atmospheric neutrino 
anomaly observed earlier by the Superkamiokande detector j2j. These, as well as the 
solar neutrino deficit measurements reported in a number of experiments jSl IH El E] 
yield valuable information on the neutrino mass differences and mixing angles [7j. 
The simplest scheme which accounts for these results is that in which there are 
just three light neutrino mass eigenstates with a mass hierarchy analogous to the 
quarks and charged leptons, and the observed phenomena of neutrino oscillations 
can be accommodated by a mixing matrix in the lepton sector — analogous 
to the well-studied quark rotation matrix [H]. If, however, the LSND result ^| is 
confirmed, it would imply a fourth, sterile, neutrino , separated in mass from the 
other neutrinos by typically 0.4-=- 1 eV. In that case, there might be even more such 
(sterile) neutrinos. 

While impressive, and providing so far the only evidence of new physics, the 
solar and atmospheric neutrino experiments do not probe the nature of the neutrino 
masses, i.e., they can not distinguish between the Dirac and Majorana character of 
the neutrinos. The nature of neutrino mass is one of the main unsolved problems 
in particle physics and there are practically no experimental clues on this issue. If 
neutrino are Dirac particles, then their masses can be generated just like the quark 
and charged lepton masses through weak SU(2)-breaking via the Yukawa couplings, 
T^D = hyv/ a/2, where hy is the Yukawa coupling and v = a/2(0°) = 246 GeV is the 
usual Higgs vacuum expectation value. In that case, to get < 1 eV neutrino masses, 
one has h^, < 10^^^, which raises the question about the extreme smallness of the 
neutrino Yukawa coupling. 

If neutrinos are Majorana particles then their mass term violates lepton number 
by two units AL = ±2 ^jjj. Being a transition between a neutrino and an antineu- 
trino, it can be viewed equivalently as the annihilation or creation of two neutrinos. 
In terms of Feynman diagrams, this involves the emission (and absorption) of two 
like-sign l^-boson pairs {W~W~ or W^+W^+). If present, it can lead to a large num- 
ber of processes violating lepton number by two units, of which neutrinoless double 
beta decay (PPou) is a particular example. The seesaw models 12j provide a natural 
framework for generating a small Majorana neutrino mass which is induced by mix- 
ing between an active (light) neutrino and a very heavy Majorana sterile neutrino 
of mass Mat. The light state has a naturally small mass m^, ~ mj^/MN <^ rriD, 
where is a quark or charged lepton mass. There is a heavy Majorana state 
corresponding to each light (active) neutrino state. Typical scale for in Grand 
unified theories (GUT) is of order the GUT-scale, though in general, there exists a 
large number of seesaw models in which both rriD and M^v vary over many orders 
of magnitude, with the latter ranging somewhere between the TeV scale and the 
GUT-scale 

If Mtv is of order GUT-scale, then it is obvious that there are essentially no 
low energy effects induced by such a heavy Majorana neutrino state. However, if 
Mtv is allowed to be much lower, or if the light (active) neutrinos are Majorana 
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particles, then the induced effects of such Majorana neutrinos can be searched for in 
a number of rare processes. Among them neutrinoless double beta decay, like-sign 
dilepton states produced in rare meson decays and in hadron-hadron, lepton-hadron, 
and lepton-lepton collisions, and e — *• /^-conversions have been extensively studied. 



(See, e.g., the papers: /?/5o^ Hi EE!, K+ ^ n-fi+i2+ [HI UHl Cni EDI EH El, 

pp e^e^X pp ^ i^i^W^X 121, e±p ^ uU^i'^X [23 EE], the nuclear 



/i~ — i> e+ 123 ESj and fi^ — > /i"*" [29^ conversion.) 

Of the current experiments which are sensitive to the Majorana nature of neu- 
trino, the neutrinoless double beta decay, which yields an upper limit on the ee 
element of the Majorana mass matrix, is already quite stringent. The present best 
limit posted by the Heidelberg-Moscow experiment [Hj is: {rriee) = \ J2iViUei^i\ < 
0.26 (0.34) eV at 68% (95%) C.L., where rji is the charge conjugation phase factor of 
the Majorana neutrino mass eigenstate Ui = rjivf. Despite some dependence of the 
actual limit on the nuclear matrix elements, this limit severely compromises the sen- 
sitivity of future e~e~ colliders |SIlE2l5 as well as of searches in the e"e~ final states, 
such as pp e~e~X, induced by a Majorana neutrino, discussed here. Hence, it 
is exceedingly important to push the /5/?oi^-frontier; currently there are several pro- 
posals being discussed in literature, which will increase the sensitivity to {niee) by 
one to two orders of magnitude [13 CHI EH, with the GENIUS proposal reaching 
(mee) ~ 10~^ eV HS|. Likewise, precision electroweak physics experiments severely 
constrain the mixing elements, namely J2n \UeN\'^, with £ = e, /i, r jHHl Ell EH] ■ 

Taking into account these constraints, we investigate in this paper the sensitivity 
to the Majorana neutrino induced effects involving same-sign dilepton production. 
We work out the following two processes in detail: 

(i) Dilepton production in the high-energy proton-proton collision 



for M = K, D, Ds, B. 

We obtain discovery limits for heavy Majorana neutrinos involved in the process 
at the LHC energy ^/s = 14 TeV. Using the present limits on the branching 

ratios of rare decays (j21) we set the upper bounds on effective Majorana masses. 

From the existing bounds on the elements of the effective Majorana mass matrix, 

the indirect constraints on the branching ratios in question are deduced. 

2 Dilepton production in high-energy pp collisions 

We have calculated the cross section for the process (HJ at high energies. 



pp ti'-^x 




with ^, £' = e, fi, T at LHC; and 

(ii) in rare meson decays of the type 



M+ 



(2) 




(3) 
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via an intermediate heavy Majorana neutrino N in the leading effective vector- 
boson approximation [32] neglecting transverse polarizations of W bosons and quark 
mixing. We use the simple scenario for neutrino mass spectrum 

and single out the contribution of the lightest Majorana neutrino assuming 

a/s <t: 171^2- 

The cross section for the process in question is then parameterized by the mass m^v 
and the corresponding neutrino mixing parameters U£n and Ue'N'- 



a (pp ^ te+x) =C[l- -5u'j \UM>N? F {E, m^) , (4) 

with 



and 



F{E,mN) = { / — / — / — p{x,xs)p [-,-s] h \ - ] w { 

\mw/ Jzo z Jz y Jy X \x X / \y ) \ 

Here, Zq = Am^r/s, E = y/s, and 



m 



(5) 



N 



w(t) = 2 + \ In (t + 1 

^' t + 1 t{t + 2) ^ ' 

is the normalized cross section for the subprocess W^W^ t^^^ (in the limit 
it is obtained from the well-known cross section for e~e~ — > W~W~ |31j using 
crossing symmetry). The function h{r) defined as 

/i(r) = - (1 + r)lnr - 2 (1 - r) 

is the normalized luminosity (multiplied by r) of W~^W~^ pairs in the two-quark 
system jHE!, and 



p (x, Q^) = x^Qi (x, Q^) = x(u + c + t + d + s + bj 

i 

is the corresponding quark distribution in the proton. 

In the numerical calculation of the cross section the CTEQ6 Fortran codes 
for the parton distributions [H] have been used. The reduced cross section (0) 
as a function of the neutrino mass m-Tv is shown in Fig. ^ for the LHC energy 
^fs = lA TeV. 

We assume the mixing constraints obtained from the precision electroweak data 

m 

E \UeN\^ < 6.6 X 10'3, E \U^,n\^ < 6.0 X 10-3 (1.8 x lO'^) , 

E|f/rjv|^ < 1.8 X 10-2(9.6 X 10-3). (6) 
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Figure 1: Left: The reduced cross section F{E,mN) defined in the text for dilepton 
production as a function of the heavy Majorana mass at LHC with E = 14 TeV. 
Right The same as the left figure but for hghter Majorana neutrinos. 



The bound on the mixing matrix elements involving fermions depends on the un- 
derlying theoretical scenario. A mixture of known fermions with new heavy states 
(here, a Majorana neutrino) can in general induce both flavour changing (FC) and 
non-universal flavour diagonal (FD) vertices among the light states. The FC cou- 
plings are severely constrained for most of the charged fermions by the limits on 
rare processes |HH]- In Refs. [HUES], the FD vertices are constrained by the elec- 
troweak precision data, which we shall use here. There are two limits obtained on 
the FD vertices, called in jHHinS] the single limit and joint limit, obtained by allow- 
ing just one fermion mixing to be present or allowing simultaneous presence of all 
types of fermion mixings, respectively. The resulting constraints are more stringent 
in the single limit case (see numbers in parentheses on the r.h.s. in Eq. (^) than 
in the joint limit case, where the constraints are generally relaxed due to possible 
accidental cancellations among different mixings. In our analysis, we shall use the 
conservative constraints for the joint limit case. 

We must also include the constraint from the double beta decay PPou, mentioned 
above. For heavy neutrinos, ttin ^ 1 GeV, the bound is [SI] 

< 5 X 10"^ TeV-\ (7) 

In calculating the cross sections for the ir and rr processes, we have used the 



2 Vn 



N{heavy) 
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effective value 

\UrN\lff = B.^ \U^Nf < 3.1 X 10-3 (8) 

with B,-^ = Br (r" — » ii^V^Ut) = 0.1737 [IH], as this r-decay mode is most suitable 
for the like-sign dilepton detection at LHC (see, e.g., |26j). 

Combining the constraints of Eqs. (jHI), and ((H) and taking a nominal lu- 
minosity L = 100 fb"^, we obtain aL < 1 for all ii' channels. Therefore at 
the LHC there is no room for observable signals for same-sign dilepton processes 
pp —>■ ii'X [i, i' = e, /i, r) due to the existing constraints for the mixing elements 
\Uei\f\ from the precision electroweak data and neutrinoless double beta decay. 

Let us take into account a recent proposal to increase the instantaneous LHC 
luminosity £ to a value of 10'^^ cm-^s"^ [SHI, i.e., a total luminosity L = Cxi year ~ 
3200 fb-\ 

Using the upgraded LHC luminosity and demanding n = 1, 3 events for discovery 
(i.e., aL > n), we present the two-dimensional plot for the discovery limits in Fig. |21 
for the case of identical same-sign leptons [i = i'). 




^ n=l ^ e X PPov |i I 

* n=l ■ n=3 ^ e -»- T -»- [Sppv |i 



Figure 2: Left: Discovery limits for pp — > i'^i'^X as functions of rriN and |f/w|^ for 
y/s = U TeV, L = 3200 fb^^ and various values of n, the number of events. We 
also superimpose the experimental limit from /S/So^-decay (Eq. ((Tj)), as well as the 
experimental limits on |f/£Ar|^ [horizontal lines for i = e, n (Eq. ® ), and r (Eq. ©)]. 
Right The same as the left figure but for lighter Majorana neutrinos. 

Discovery limits for the case of distinct same-sign leptons, ii' = efi, er, fir, are 
shown in Fig. 01 
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Figure 3: Left: Discovery limits for pp — > i^£'~^X, ii' = e/i, er, fir as functions of 
rriN and \/2 \U£NUi/]y\ for y/s = 14 TeV, L = 3200 fb^^ and various values of n, the 
number of events. We also superimpose the limits on ^/2\UeNUi'N\ obtained from 
the experimental limits [Eqs. ((Tj), (0), and (jH))]. Right The same as the left figure 
but for lighter Majorana neutrinos. 



/^From Figs. |2 and 121 we see that the existing strong constraints for the mixing 



2 

elements \U£n\ allow a possibility to observe only the same-sign /i/i and /ir processes 
after increasing the nominal LHC luminosity L = 100 fb~^ by a factor of about 30. 
For this case, LHC experiments will have a sensitivity to heavy Majorana neutrinos 
of mass niN ^ 1-5 TeV. 

Before concluding this section, we would like to add a comment on Ref. 
where ostensibly new and improved bounds on the effective Majorana masses 



N 



(9) 



have been obtained using the HERA data on ep collisions. The authors in Ref. [2^] 

have assumed that the cross sections for the processes e^p — W i^i'^X are pro- 
portional to {mu')'^. This, however, is true only for light Majorana neutrinos. For 

the heavy Majorana neutrino case, the cross sections depend on the factor (m^J) , 
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where 



TJln 



(10) 



i.e. the effective inverse Majorana masses. As a consequence, Eq. (6) of Ref. [SE] 
does not give new physical bounds on the hght Majorana neutrino masses (as the 
resulting cross section is too small), and for the heavy Majorana neutrino case, their 
formula is not applicable (see also comments on Tables El and El below) . Hence, 
contrary to the claims in Ref. [SHl , HERA data does not place any new limit on the 
Majorana mass matrix. 



3 Rare meson decays M+ M'-i+i'+ 

We now take up rare meson decays of mesons of the type (j2I) mediated by Majorana 
neutrinos. We shall take the mesons in the initial and final state to be pseudoscalar. 
The lowest order amplitude of the process is given by the sum of the tree and the box 
diagrams shown in Fig.j^taken from Ref. [221 (for earlier work, see Refs. |221EICS1)- 




Figure 4: Feynman diagrams for the rare meson decay — M'~i^i'^. Here 
is a Majorana neutrino, bold vertices correspond to Bethe-Salpeter amplitudes for 
mesons as bound states of a quark and an antiquark. There are also two crossed 
diagrams with interchanged lepton lines. 

It is well-known that the tree diagram amplitude can be expressed in a model 
independent way in terms of the measured decay constants of the pseudoscalar 
mesons in the initial and final state, /m and J'm'- On the other hand, the box diagram 
depends in general on the details of hadron dynamics. In the usual folklore, the tree 
diagram dominates and is often used to set limits on effective Majorana masses 
pT] . Calculations based on various quark models for mesons have shown that this 
dominance really holds in the decay n^fi^fi^ for rather small neutrino masses 

PTj . It can be explained, at least partly, by a colour suppression factor 1/Nc = 1/3 
present in the box amplitude. Indeed, we see in Fig. that in the quark loops in the 
tree diagram the colour summation takes place independently in the two loops. Not 
so in the box diagram, where the colour of the quark qi (antiquark gg) must be the 
same as the colour of the quark (antiquark q^) in the lower (upper) part of the 
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box diagram. We also note that in some cases there is Cabibbo suppression of the 
tree or box amphtude due to smallness of the corresponding CKM matrix elements. 

In this paper, we calculate branching ratios for rare meson decays Q for two 
limiting cases of heavy and light Majorana neutrinos. In particular, we find that for 
the case of heavy neutrinos, rriN ^ rnM, the box contribution to the decay amplitude 
can also be expressed through /m and /m' independent of a specific structure of the 
Bethe-Salpeter vertex for the meson in question. 

The width of the rare decay M+(P) M'~{P')i+{p)i'+{p') is given by 

T,,= (l--5,,)ji2.) 6^HP+P + P-P) 23 (2.) V - (1^) 

Here At (A^) is the tree(box)-diagram amplitude expressed in the Bethe-Salpeter 
formalism of Ref. jlH] as 

^^ = 7^8/ d^d^'H^m^ , (12) 



(27r 

where the lepton tensor is given by 

4 n^"* n^l^ 



4 pf - mlr pf - 



X v'^ (p) 



lalf3 , 7/37. 



N 

1+7^ 



{Pi -P) -^N (Pi - P') - ^ 



N 



v{p') , (13) 



and r/Tv are the charge conjugation phase factors; i = t, b, implying tree and box 
contributions, respectively, and 

Pt = P, p[ = P'; Pt=liP-P') + q'-q, p', = \{P - P') - q' + q. 
The hadron tensors 

H^l = Tr (g) ^127.^} Tr (g') ^437.^^} , 

<^ = Tr (g) ^137.^^ {q') ^42%^} (14) 

are expressed in terms of the elements of the CKM matrix, Vjk (the subscripts on 
V correspond to the quark line numbering in Fig. 0]), and the mo del- dependent 
Bethe-Salpeter amplitudes for the mesons 



(g)=7ni-5M/^)^(P,g) <^G, 

where 6m = (^i + m2)/M'^, M is the mass of the meson having a quark qi and 
an antiquark ^1,2 are the quark masses, q = {pi — P2)/2 is the quark-antiquark 
relative 4-momentum, P = pi + p2 is the total 4-momentum of the meson; the 
function ip{P, q) is model dependent, and 0g is the SU (Nf) x S'?7(Ai'c)-group factor. 
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Table 1: The pseudoscalar decay constants Jm for the indicated mesons 



Meson 


/m [MeV] 




T ''7 

130.7 




159.8 


D+ 


228 


Dt 


251 


B+ 


200 



The decay constant /a/ of the meson M is expressed through the amphtude ^is 
folfows: 

ifuP' = (O \q, (0) (0)| M (P)) = -^^iVc / ^Tr {7'^7 V (g)} , 

/M = 4yAr:5A/ / ^^(P,g), (15) 

where the sum over cofour indices is imphed. The values of the pseudoscalar coupling 
constants /m experimentally measured (for vr and K jHH]) and the central values 
calculated using lattice QCD (for D, Ds, and B mesons ^T]) are shown in Tabled 
We have neglected the errors on these quantities, as we shall see that this will 
not compromise our conclusions in any significant way. For all mesons in question, 
niM ^ T^w^ and we can use the leading current- current approximation in the lepton 
tensors (fT3|) . 

Below we consider two limiting cases of heavy and light Majorana neutrinos. 



Heavy neutrinos: itln ^ tum 

For this case, the tree and box lepton tensors are equal to each other in the leading 
order of the expansion in 1/m^: 

Lr = ^' = -16g^''L{p,p') , 

L{p,p') = G'^Y^UmUeNVN— (^{p) ^^^p')) • (16) 

/^From Eqs. (fT^. (Hlj), (US)), and (fT^ we obtain the total amplitude of the decay 

A = At + A, = AKvfMfM'iP ■ P')L{p,p') , 

Kv = VuV^s + -^VM2 , (17) 

which is model independent in this limit. 
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Using Eqs. (fTTjl and (fTTjl we calculate the decay width: 



Here the effective inverse Majorana neutrino mass (rn^^j is introduced in the form 
of Eq. (fTn|) and is the reduced phase space integral. For identical leptons 



21 



4;zn^ 

= / dz{z- 2zo) ' ' 



(l - ^) - .) {z, - z) 



\^z^-zf. (19) 

4zo 



For the case of I and ^' being distinct leptons, assuming vriii jmi <^ 1, in the leading 
approximation, we have 

=2 T — iz- Zof [{Z^ - Z) {Z2 - z)f^ il + Z3- zf , (20) 



where the variable of integration is 2; = (P — P') / m\j and the parameters Zk are 
defined as 

2 2 

_ f^H _ /-I / — \2 _ /-I ; x2 _ ""^M' 

mil 

Using Eqs. (fTHj). (fT^ . (plj) . Table [T] and the experimental values for the total 
decay widths of mesons from ||38| , we have calculated the branching ratios 



R (An r (M+ M-f+f +) 

^'"^^^ - V (M+ all) ■ ^^^^ 

Comparing the results with experimental bounds on Bn' taken from (for K 
decays), ^ (for D and decays), and [38j (for D+ K-i+i'+, B+ 7r-£+f+, 
and — >■ K~i~^i'~^ decays), the upper bounds for the effective inverse Majorana 
masses fllOj) have been obtained. The results are shown in Table |21 

We can also obtain the indirect upper bounds on the branching ratios using the 
constraint ((Zj) from (3(3ou on the element of the effective inverse mass matrix. 

For other elements, assuming one heavy neutrino scenario, we set 

mji,}) < (84.1 GeV)-\ 



using the current mass limits on neutral heavy leptons of the Majorana type 
The corresponding indirect bounds are shown in the last column of Table |21 
Our estimate for the K decay branching ratio, 

B^^{K) =b[k+ ^ 7r-/i+/i+) = 2.5 x 10-^° MeV^ ■ {m~'^y (22) 

should be compared with the corresponding updated result of Dib et al. ■ Tran- 
scribed to our notations, the relevant decay width is given by the following expression 
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Table 2: Bounds on (m^^j and indirect bounds on the branching ratios B^^, (M) 
for the rare meson decays M+ — > M'~l'^(!^ mediated by Majorana neutrinos (with 
mjss^ rau) and present experimental bounds 



Rare decay 


bxp. upper bounds 


Ineor. estimate tor 


Bounds on 


T J 1 J 

Ind. bounds 




on B^^,(M) 


B«,(M)/(m7/) [MeV2] 


(mie') [keV] 


on B^^,(M) 




6.4 X 10-^" 


8.6 X 10-^" 


1200 


2.2 X 10-^" 




3.0 X 10-9 


2.5 X 10-1° 


300 


3.5 X 10-2" 




5.0 X 10"^" 


8.4 X 10-1" 


1300 


1.2 X 10-iy 


D+ Tr-e+e+ 


9.6 X 10"^ 


2.2 X IQ-y 


4.8 


5.5 X 10-^" 


D+ Tr-iJ+fi+ 


1.7 X 10-^ 


2.0 X IQ-'^ 


11 


2.8 X 10-19 


D+ 7r-e+/Lt+ 


5.0 X 10"^ 


4.2 X IQ-y 


9.2 


5.9 X 10-19 


D+ K-e+e+ 


1.2 X 10"^ 


2.2 X IQ-y 


4.3 


5.5 X 10--^" 


D+ K-^i+^i+ 


1.2 X 10"^ 


2.1 X IQ-y 


4.1 


3.0 X 10-19 


D+ K-e+n+ 


1.3 X 10~^ 


4.3 X IQ-'^ 


5.7 


6.1 X 10-19 


D+ 7r-e+e+ 


6.9 X 10-^ 


1.9 X 10-« 


5.2 


4.8 X 10-29 


D+ ^ 7r-fi+fi+ 


8.2 X 10"^ 


1.8 X 10-« 


15 


2.5 X 10-1^ 


D+ 7r-e+/i+ 


7.3 X 10-4 


3.6 X 10-« 


7.1 


5.1 X 10-1^ 


D+ K-e+e+ 


6.3 X 10-4 


2.2 X lO-'^ 


1.9 


5.5 X 10-^" 


Dt ^ K-^^+^i+ 


1.8 X 10-4 


2.0 X IQ-'^ 


3.4 


2.8 X 10-19 


D+ K-e+ti+ 


6.8 X 10-4 


4.2 X IQ-y 


2.5 


5.9 X 10-19 


B+ 7r-e+e+ 


3.9 X 10-^ 


(0.3^ 1.9) X lO-'^ 


0.3^0.7 


4.8 X 10-^" 


B+ 


9.1 X 10-^ 


(0.3^ 1.9) X lO-'^ 


0.2^0.5 


2.7 X 10-19 




6.4 X 10-^ 


(0.6^3.8) X lO-'^ 


0.3^0.8 


5.4 X 10-19 


D+ — --T+T- 




(0.2 ^ 1.2) X 10"'' 




1.7 X 10-19 


B+ Tr-e+T+ 




(1.0^6.2) X 10-1" 




4.0 X 10-19 


B+ TT-fl+T+ 




(1.0^6.2) X 10-1" 




8.8 X 10-2" 


B+ K-e+e+ 


3.9 X 10-^ 


(0.2^ 1.5) X 10-1" 


0.07^0.20 


3.8 X 10-^1 


B+ 


9.1 X 10-^ 


(0.2 H- 1.5) X IQ-i" 


0.05 H- 0.13 


2.1 X 10-2U 


B+ K-e+^+ 


6.4 X 10-^ 


(0.5 H- 2.9) X IQ-i" 


0.09-=- 0.21 


4.1 X 10-2" 


B+ -> K-T+T+ 




(1.3^8.4) X 10-1^ 




1.2 X 10-21 


B+ K-e+T+ 




(0.7 -=-4.4) X 10-11 




6.2 X 10-21 


B+ K-ij,+T+ 




(0.7 H- 4.4) X IQ-ii 




6.2 X 10-2i 
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yielding a branching ratio 

We note that including the box diagram gives the correction factor (see Eq. (fTTjl ) 
(1 + l/Nc)"^ = 16/9, i.e. the numerical coefficient 1.6 in the above equation gets 
replaced by 2.8, which is close to the numerical coefficient 2.5 in our Eq. ()22j) . 
In addition, a rough estimate of the branching ratio, 

obtained in ^Hj (and conffimed in ^2Qj) is also in agreement with Eq. (j22I)- 

^From Table El we see that the present experimental bounds on the branching 
ratios of the rare meson decays are too weak to give interesting bounds on the 
Majorana mass. From them and Eq. ()18|) . we see that the direct bounds on the 

effective masses (rnj^'^ of heavy Majorana neutrinos are much sma//er than the 
difference of meson masses, tum — ■ The indirect bounds estimated on B«/, as 
discussed above, and given in the fifth column are so small that they can not be 
realistically tested in any current or planned experiments. 



Light neutrinos: trn ^ fng , 



In this case, assuming the tree diagram dominance, the lepton tensor of Eq. ()1H|1 . 
neglecting m% in the dominators, can be approximated as 



N 



+ 



(Pk - p) {pk - P') 



i + r 



v{p') 



Using Eq. (fTT|) with A At, we obtain the width of the decay (0) for the case of 
light Majorana neutrinos 



GpTTlj^yj 2 n2 IT/ T/ |2 \2 J 



167r^ 



(23) 



where {mui) is the effective Majorana mass (jH)). Here the phase space integral 0^/ 
has a rather complicated expression but in the realistic limit of massless leptons, 
mi/niM — > and mii /itlm — 0, it can be approximated as 



1 - \^u'^ V (^3) ; V (^3) = dzz [{zi - z) {z2 - z)] 

2^3 + ^ (1 - ^3)' 
6 



1/2 



;i - ^3) 



+ 2:3 (1 + 2:3) In 2:3, 



where Zk are the same as in Eq. ((201). 

Using Eq. ((221); "^^ have calculated the branching ratios (PT|) and obtained the 
direct upper bounds on the elements of the effective Majorana mass matrix. The 
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indirect limits on the branching ratios have been also obtained with use of a rather 
stringent constraint on the ee element, 

(mee) < 1.0 eV, 

liiom the jSPoi, Heidelberg-Moscow experiment ^3] (see also comments in and 
a weaker constraint on other matrix elements used in [221 

{mu>) < 9 eV , 

which has been deduced in the three light neutrino scenario assuming the upper 
bound of 3 eV on the mass of the known neutrino [IHl (more stringent but model 
dependent bounds on (m«,) have been obtained in |Tfj). 

Our results are shown in Table 01 We note that for the K~i~^i'~^ decays, 

the tree diagram is strongly Cabibbo suppressed and the box diagram must be 
included even for the case of light neutrinos. We have obtained rough estimates 
of these decay widths assuming that the reduced (i.e., without CKM and colour 
factors) tree and box amplitudes are of the same order and replacing in Eq. 
the factor iVuVis] hj \Kv\ (see Eq. (HH)). It gives a numerical correction factor 
of about 57. 

We conclude once again that the present experimental bounds on the branching 
ratios of the rare meson decays are too weak to set reasonable limits on the effective 
masses of light Majorana neutrinos. 

As for the heavy Majorana neutrino case, one of our results, 

B ^ 7r"/i+/i+) = 1.4 X 10"^° MeV-^ ■ {m^^f , (24) 
is close to the corresponding result (in our notations) of Ref. |221, 

y(DGKS) ^ ^-^+^+) = 4.0 X 10-31 . (m^^)' , 

or 

B^^^n^) = 1.5 X 10-° MeY'.{m-l)\ 
and the minimal value of a rough estimate of Refs. [T?H I2()j. 

B^^(K) ~ 0.2 X 10-(i3±2) (^^^^^ /loO MeV)' , 

gives the same order of magnitude of the branching ratio as Eq. ()24j) . 

^From Eq. we obtain the constraint (see Table El) {fn^i^j.) < 470 GeV, which 
is almost the same as the one obtained in |2I1: {m^^) < 500 GeV. But we have to 
stress, in contrast to the conclusion of Ref. [23] , that there is no reasonable limit at 
all that emerges /^from rare decays, as Eq. (f^^ is valid for m^v <^ fn^ ~ 100 MeV, 
and since \U^n\ < 1 by unitarity, the obvious inequahty (see Eq. Q) {m^^) < 
\J2n''tT'n\ holds. We note that Ref. [2] was also criticized on the same grounds in 

Refs. [201 1221- 
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Table 3: Bounds on {mui) and indirect bounds on the branching ratios B^^, (M) 
for the rare meson decays — *■ M'~£+£'"'" mediated by Majorana neutrinos (with 
m-Tv <^ m£,m^/) and present experimental bounds 



Rare decay 


Exp. upper bounds 
on B„,(Af) 


Theor. estimate for 
B^^,(M)/(m«,)2 [MeV-2] 


Bounds on 
(m«,) [TeV] 


Ind. bounds 
on B„,(M) 




TT! 

6.4 X 10"^'^ 


rrn 

5.1 X IQ-^" 


0.11 


fTTi 

5.1 X 10-^^2 




3.0 X 


1.4 X IQ-^u 


0.47 


1.1 X lO-^'^ 




T7T 

5.0 X 10"^" 


^TTi 

6.2 X 10-^^ 


0.09 


^TTi 

5.0 X IQ-^'^ 


D+ Tr-e+e+ 


y 

9.6 X 10"^ 


m 

1.2 X 10-^^ 


280 


n — 

1.2 X 10-^^ 


D+ 7r-/i+^+ 


p 

1.7 X 10"^ 


n 

1.2 X 10-^^ 


120 




9.7 X 10-^^2 


D+ 7r-e+/i+ 


y 

5.0 X 10"^ 


m 

2.3 X IQ-^^ 


150 


n 

1.9 X 10-^^ 


D+ K-e+e+ 


1.2 X 10"^ 


2.3 X 10-21 


230 


2.3 X 10-33 


D+ K-^i+n+ 


1 

1.2 X 10-4 


rn 

2.3 X 10-2^ 


230 


TT 

1.8 X 10-31 


D+ K-e+n+ 


1.3 X 10-4 


4.5 X 10-2^ 


170 


3.7 X 10-31 


D+ TT~e+e+ 


6.9 X 10-4 


1.5 X IQ-'^^ 


210 


1.5 X 10-3^ 


Dt ^ 


y 

8.2 X 10-^ 


rrri 

1.5 X 10-^^ 


74 


TTTi 

1.2 X 10-3*^ 


Df Tr-e+ iJ,+ 


1 

7.3 X 10-4 


TTTi 

3.1 X IQ-^^ 


150 


TTTi 

2.5 X 10-3*^ 


D+ K-e+e+ 


1 

6.3 X 10-4 


TTiT 

5.6 X 10-22 


1100 


m 

5.6 X 10-34 


D+ ^ K-^^+^,+ 


1 

1.8 X 10-4 


TTTT 

5.6 X 10-22 


570 


4.5 X 10-32 


D+ K-e+^+ 


1 

6.8 X 10-4 


rn 

1.1 X 10-21 


780 


8.9 X 10-32 


B+ -> TT-e+e+ 


3.9 X 10--^ 


(0.3^ 1.8) X 10-2^ 


(1.5- 


- 3.6) X 104 


1.8 X 10-35 


B+ Tr~ ^1+11+ 


9.1 X 10--^ 


(0.3^ 1.8) X 10-2^ 


(2.2- 


- 5.5) X 104 


1.5 X 10-3^ 


B+ TT-e+n+ 


6.4 X 10-^ 


(0.6^3.6) X 10-2^ 


(1.3- 


- 3.3) X 104 


2.9 X 10-3^ 


B+ -K-T+T^ 




(1.5^9.6) X 10-25 




7.8 X 10-35 


B+ 7r-e+r+ 




^ ; ?n 

(0.4^ 2.4) X 10-2^^ 




n — 

1.9 X 10-33 


B+ 7r-//+r+ 




^ ; ?n 

(0.4^ 2.4) X 10-2^^ 




n — 

1.9 X 10-33 


B+ K-e+e+ 


n 

3.9 X 10-'^ 


(0.2^ 1.2) X 10-24 


(0.6^ 1.4) X 10^ 


— 

1.2 X 10-35 


B+ K-fi+fi+ 


9.1 X 10-'^ 


(0.2^ 1.2) X 10-24 


(0.9^2.2) X 10^ 


9.7 X 10-35 


B+ K-e+fi+ 


6.4 X 10--^ 


(0.4^ 2.4) X 10-24 


(0.5^ 1.3) X 10^ 


1.9 X 10-34 


B+ K-T+T+ 




(1.0^6.1) X 10-25 




4.9 X 10-35 


B+ K-e+T+ 




(0.2^ 1.2) X 10-24 




9.7 X 10-35 


B+ K-H+T+ 




(0.2^ 1.2) X 10-24 




9.7 X 10-35 



4 Conclusion 

We have examined two processes mediated by Majorana neutrinos: the production 
of like-sign dileptons £~^£'~^ in proton-proton collisions at the LHC energy and in 
the rare decays of K~^, D^, Df, and mesons of the type —>■ M'~[^£'^ . 
We find that at LHC for the nominal luminosity L = 100 fb~^ there is no room 
for observable same-sign dilepton signals due to the existing constraints from the 
precision electroweak data and neutrinoless double beta decay. But increasing the 
nominal LHC luminosity by a factor of about 30 will allow a possibility to observe 
the same-sign /i/i and /ir processes mediated by a heavy Majorana neutrino of mass 
niN ^ 1.5 TeV. Data from HERA do not have an impact on the Majorana mass 
matrix — despite claims to the contrary jSEj- However, precision electroweak data 
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may lead to more severe constraints on the fermion mixing angles than worked out 
in EH] and used by us. As for the rare meson decays, present direct bounds on 
their branching ratios are too weak to set reasonable limits on effective Majorana 

masses, (m«') (for light neutrinos) and (mj^l'j (for heavy ones). Therefore, to 
have an impact on the Majorana mass matrix, a very substantial improvement of 
the experimental reach on the lepton-number violating rare meson decays is needed. 
Conversely, if a same sign dilepton signal is seen in any of the meson decay channels 
listed in Tables El and El in foreseeable future, it will be due to new physics other than 
the one induced by Majorana neutrinos, such as R-parity violating supersymmetry 
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